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Abstract 

Background: Two important parametrizations of momentum-dependent nucleonic fields, pro- 
posed for the simulations of central heavy-ion collisions, one by Gale et al. and the other by Welke 
et al., suffer from practical limitations. The first gives rise to mean fields isotropic in momentum, 
even when underlying momentum distributions are anisotropic, making descriptions of early non- 
equilibrium stages of collisions unrealistic. The second parametrization gives rise to anisotropic 
mean fields, but is computationally expensive, since the mean field has to be computed separately 
for every location of a nucleon in phase-space, through folding. 

Purpose: Here, we construct a parametrization of the nucleonic mean field, that yields an 
anisotropic mean field for an anisotropic momentum distribution and is inexpensive computa- 
tionally. To demonstrate the versatility of our parametrization, we take the case of results from 
CN . the parametrization by Welke et al. and attempt to approximate them. 

, Method: In arriving at a suitable anisotropic mean field potential, we draw, on one hand, from 

the idea behind the parametrization of Gale et al, of a separable expansion of the potential energy, 
and, on the other, from the idea of a parallel expansion of the energy and mean field, in anisotropy. 
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Results: We show that using our novel parametrization we can qualitatively and partially quan- 
titatively reproduce the features of the mean-field parametrization of Welke et al. . 



Conclusions: This opens up the possibility of exploring the effects of mean-field anisotropy in 
collisions, without the penalty of computational cost behind the folding parametrization. 
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I. INTRODUCTION 



In heavy-ion collision theory, the interplay of the momentum dependence of the nucleon 
optical potential U and collision observables is of utmost importance. Not only does mo- 
mentum dependence play a significant role in generation of collective flow, according to the 
transport calculations [1], but it is also crucial for particle production [2]. In particular, in 
order to properly constrain the nuclear compression modulus K, defined as 



where Pf is the Fermi momentum and E/A the binding energy per nucleon, one must take 
the momentum dependence of nucleon-nucleon interactions into account. In 1976, Blaizot et 
al. [3] showed that K could be inferred by measuring the energy of the isoscalar monopole 
resonance in medium and heavy nuclei, and arrived at the value of K — 210 ±30 MeV. Flow 
data, however, for quite some time seemed to be describable by both a momentum-dependent 
"soft" eos (K « 210 MeV) or a momentum-independent "hard" eos (K « 380 MeV) [4, 5]. 
Eventually, though, it became possible to decide on the momentum dependence with heavy- 
ion observables. Pan and Danielewicz [1] demonstrated that the best agreement between 
sideward flow data and the results from transport calculations could be obtained when 
applying a momentum-dependent "soft" eos. 

Several approaches to the momentum dependence of the optical potential have been 
put forward in the literature (see [6] for a review article). The particularly well-known 
parametrizations, utilized directly or indirectly in Boltzmann-Uehling-Uhlenbeck (BUU) [7] 
calculations, are the early parametrization by Gogny [8], the ansatz by Gale, Bertsch and 
Das Gupta [5] and the one by Welke et al. [9]. In Sec. II, the latter two parametrizations 
are discussed in more detail, and the results by Welke et al. will serve as a reference for our 
calculations. Following Landau quasiparticle theory [10], inherent in these approaches are 
a functional expression for the potential energy density V and for its functional derivative 
with respect to the single-particle phase-space density /(r, p), 
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namely the nucleonic mean field (or optical potential) U. 

If the density / is anisotropic in momentum space, one would expect, by symmetry 
considerations, the mean field to exhibit anisotropy as well. Parametrizing / by some 
anisotropy scale e can be useful in discussing partially equilibrated scenarios in momentum 
space, that develop after the nucleons from the original target and projectile Fermi spheres 
had a chance to collide with the nucleons from the counterpart nucleus. Thus, the idea of 
making the mean field U itself anisotropic for the scenario of anisotropic / has great appeal 
to us, especially since some of the models employed in practice lack this feature. 

Our work (based on the M.Sc. Thesis [11]) is devoted to deriving a parametrization for the 
energy density V and the optical potential U in heavy-ion collisions for anisotropic systems, 
that can be efficiently used in BUU calculations, such as within the code by Danielewicz 
[12], exploring different anisotropies of U. The construction involves introducing separable 
interactions in p-space where different terms represent, in particular, different spherical 
harmonics through which anisotropy explicitly enters the mean field. We attempt to make 
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our approach suitable for different stages of a heavy-ion collision. Another ansatz [13-15], 
aiming at a simplified treatment of momentum anisotropies, focusses on the early stages 
of a collision, representing the momentum distributions in terms of two zero-temperature 
Fermi distributions. Since the approach by Welke et ai, in a BUU environment, requires the 
determination of a 3-dimensional integral for every relevant phase-space location at every 
time step, it is in general computationally very costly. In Sec. II, we present the models 
by Gale, Bertsch, Das Gupta (GBD) and by Welke et al. (WPKDG) and explain our own 
ansatz for V and U, based on a spherical harmonics expansion (SHE). Our SHE results 
for ground-state, excited-state and collision scenarios are represented in Sec. Ill, and our 
approach is tested there against the WPKDG parametrization. The benefits of SHE for 
BUU calculations, in terms of increased computational efficiency and reduced Monte-Carlo 
noise, which are disadvantages of WPKDG, are outlined in Sec. IV. 



II. THEORETICAL CONSIDERATIONS 



A. Implicit anisotropy in mean fields 



The special feature of GBD is that it is formulated in the local frame where the average 
nucleonic momentum (p) vanishes. Alternatively, the nucleonic momenta may be specified 
relative to (p) in another frame. This may be contrasted with WPKDG. Both models are 
based on Skyrme-type interactions [16] for the density-dependent part of the energy density 
and the mean field. In the case of GBD, the potential energy density can be written as [5] 
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Accordingly, following (2), one obtains the nucleonic mean field 
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Note that / - where not otherwise stated - represents the single-particle phase-space density, 
normalized according to p(r) = f d 3 p' /(r, p'). 

Within WPKDG, the energy density is parametrized as [9] 
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This leads to the mean field of the form 



f/wPKD G (p(r),p)= A[^-)+B f 



Po J V Po 

Inherent in both models are the 5 parameters A, B, a, C and A. From those parameters, a is 
particularly strongly tied to the incompressibility K, while A determines the high-p behavior 
of U. To constrain the parameters, a demand is placed that the energy per nucleon minimizes 
with E/A = — 16MeV at the saturation density po = 0.16 fm~ 3 . For both parametrizations, 
GBD and WPKDG, the incompressibility is chosen equal to K ~ 215 MeV, and effective 
mass ratio at saturation momentum is m*/m ~ 0.7. For WPKDG, specific values of the 
potential at p are U(p ,p = 0) ~ -75 MeV and U( y p ,p 2 /2m ~ 300 MeV) = 0. The 
GBD parametrization may be considered an approximation to WPKDG, and this issue has 
been explored to a degree in [9]. With SHE we want to develop a qualitatively better 
approximation to WPKDG, while still keeping the method computationally inexpensive, as 
GBD. 



B. Explicit anisotropy in mean fields 



The essential idea of SHE, borrowed from GBD [5], is that of a separable representation 
for the energy functional. A separable representation, with a limited number of terms, 
could be adjusted to represent results from microscopic theory, obtained for different non- 
equilibrium situations, such as Dirac-Brueckner [17]. In lieu of a microscopic theory, we 
take the intuitively appealing WPKDG parametrization and examine to what extent we 
can reproduce its results with SHE. If we can be successful here, we may be successful 
with a microscopic theory as well. Otherwise, following a phenomenological strategy, as 
common in heavy-ion collisions, we may use SHE as the starting point and see whether 
collision observables can constrain the energy functional constructed within SHE. Besides 
the separable part of the energy, generating momentum dependence of U, we employ a p- 
dependent part of the Skyrme form, such as in GBD or WPKDG. When aiming at a minimal 
number of terms in the separable expansion, while retaining the capability of describing 
mean fields anisotropic in momentum, our energy functional in SHE consists of Skyrme and 
separable scalar and tensorial quadrupole terms: 
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We do not include a dipole separable term, which we expect to be small in the center of mass. 
A dipole term is accounted for in the relativistic approach of Reference [18]. To demonstrate 
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the robustness of such a representation we will attempt to show that the results from the 
folding model WPKDG can be well approximated within our separable model, including the 
reproduction of anisotropies of U. For simplicity, we take both the scalar and the tensor term 
as symmetric in two single-particle quantities, presuming dominance of 2-body interactions. 
Inclusion of a tensor term makes our parametrization different from GBD, scalar in our 
terminology. However, we take the form of the single-particle factor to be the same as in 
GBD [5]. In the above, T a P is a second-rank Cartesian tensor [19, 20]: 
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An important feature is its tracelessness, 
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For the functional form of c(p) in (8) we choose 
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Note the parameters A; so in the scalar part of (7) and A an i SO in the tensorial part, which we 
introduced in addition to the ones inherited from GBD and WPKDG. 

According to (2), the optical potential is found by taking the functional derivative of (7) 
with respect to /. Due to the directionality involved and approximate axial symmetry of the 
phase-space density in heavy-ion collisions, the off-diagonal elements of (8) are negligible in 
practice, reducing the summation over a, /3 to diagonal elements only. Here the tracelessness 
(9) of the tensor comes into play, allowing us to express all elements in terms of T zz when 
dependence on transverse direction is weak. Upon taking a functional derivative of (7) and 
transforming to spherical coordinates we get 
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Here we used the definition of the spherical harmonic of degree 2 and order 0, 
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Therewith, the SHE optical potential (11) depends on anisotropy explicitly. The decisive 
last term vanishes for isotropic / since T zz is zero for such cases by construction. 
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TABLE I. Parameters in GBD, WPKDG and SHE models. 



Model 


A [MeV] 


B [MeV] 


C [MeV] 


a 


A 


GBD 


-144.9 


203.3 


-75.0 


7/6 


1.5^ 0) 


WPKDG 


-110.44 


140.9 


-64.95 


1.24 


1.58p<°> 


SHE 


-110.63 


132.17 


-53.04 


1.27 





Aj S o A an ; so 



1.98 1.41 



ILL RESULTS 

A. Ground-state properties 

As in GBD and WPKDG, we demand that E/A for SHE minimizes at po — 0.16 fm -3 with 
the value of E/A = —16 MeV. For consistency with the other parametrizations, we further 
demand that K ~ 215 MeV, and m*/m ~ 0.7 at the Fermi momentum of p^ = 263MeV/c 
at po- Otherwise, it would be natural to require that SHE reproduced the momentum 
dependence of WPKDG for the ground-state phase-space density 

/(r ' p)= (2W %f) ~' Pl) ' (13) 

where the degeneracy factor is g — 4. In fact, by adjusting the parameters of SHE the mo- 
mentum dependencies at po can be made nearly indistinguishable. However, then differences 
between SHE and WPKDG may be excessive in some non-equilibrium situations. Corre- 
spondingly we adopt a compromise choosing parameters so that the results have a similar 
appearance under different circumstances, but can exhibit quantitative discrepancies. 

At zero temperature, calculations of energy and optical potential can be largely done 
analytically. A couple of integrals useful for the purpose are provided in the Appendix. 
With spatial density of kinetic energy given by 

P 2 \ _ f^P' 2 



L = d »'^ p ' h (14) 



the energy per nucleon is 

E/A = J d 3 r 



I A . (15) 



£ + !7Wr),p) 



(16) 



From the energy, the incompressiblity K is calculated and, from U, the effective mass ratio 
is determined 

m* p ( d 
m m\dp 

The parameters for SHE, as well as GBD and WPKDG are provided in Table I. Figure 1 
shows the potential energy density for the three interaction parametrizations and Fig. 2 
displays the energy per nucleon. Next, Fig. 3 shows the incompressibility as a function 
of density p and Fig. 4 shows the effective mass ratio at po, as a function of momentum. 
Finally, Fig. 5 shows the optical potential as a function of momentum at the normal density. 
As we stated, by adjusting the parameters in SHE, we could obtain an excellent agreement 
between SHE and WPKDG in Figs. 4 and 5, but that at the cost of poorer agreement in 
other situations. 



6 




FIG. 1. Ground-state potential energy density as function of density of matter. 




FIG. 2. Ground-state energy per nucleon as function of density of matter. 



B. Excited nuclear matter 



To study suitability of SHE for describing situations in excited but not fully equilibrated 
matter, we consider results from SHE and WPKDG for a phase-space density which has the 
form of an anisotropic Gaussian in momentum: 
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Here, e is a parameter that regulates anisotropy of the momentum distribution that is axially 
symmetric about the longitudinal axis. Changes in e simultaneously change the longitudinal 
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FIG. 5. Ground-state optical potential as function of momentum, at normal density po- 



and transverse widths in such a manner that the density p associated with / does not change. 
Physically allowed values are e > — 1 and increasing e makes the distribution more prolate. 
One needs to relate the parameters e and cr g to other, e.g. such as temperature, in order 
to determine interesting parameter ranges. We refer to the kinetic energy density (14) to 
advance in our considerations. While an isotropic Gaussian (e = 0) yields the density 
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for an anisotropic Gaussian (e ^ 0) one obtains 
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Recalling the Boltzmann model, Gaussian standard deviation a and temperature T are 
connected via 



a' 
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(20) 



where m^ is the nucleon mass. Effective temperatures between 15 MeV < T < 170 MeV are 
of interest in heavy-ion collisions. We now require that the quantity 




l + e) 



(21) 



inferred from a comparison of (18) with (19), shall be constrained such that <r^/m N does 
not exceed the temperature bound of 170 MeV. More spread out distributions will get their 
widths quickly tempered by inelastic processes. The largest anisotropies will be permitted for 
low cr g . Taking er g equal to the width of the Fermi sphere at p = 2p , we get a g ~ 150 MeV/c 



9 



> 

CD 



-20 



-30 



-40 



^> -50 



3 

u 

a. 



-60 



-70 



-80 



(jg = 150MeV/c 
e = -0.5 



^ 0) /2 



p = 



-1 



-0.5 





COS$ 



SHE 

WPKDG 



0.5 



FIG. 6. Optical potential for an anisotropic momentum distribution (17) with e 
in the local cm., for different indicated momenta p. 



-0.5 at p = 2p , 



and arrive at —0.7 < e < 1.5. Other common sense considerations yield similar ranges of e. 
In the following we pursue the exemplary case of a g ~ 150MeV/c, on the verge of realism, 
which however allows us to consider the large range of e. The subsequent Figs. 6-10 show 
the obtained values of optical potential both in the SHE and WPKDG parametrizations, at 
different momenta p (in multiples of ground-state pp\ Sec. Ill A) in p = 2p matter with 
different momentum anisotropies, as a function of angle relative to the symmetry axis of the 
distribution. In the figures, we can see rather good correspondence between the results for 
the two parametrizations. While the SHE parametrization explicitly follows the shape of 
Y20, it is approximately also the case for WPKDG. Apparently higher multipolarity terms in 
a separable expansion for WPKDG are negligible. This is likely helped by the smoothness 
of the momentum distribution. Notably in heavy-ion collisions the momentum distributions 
tend to evolve to a smooth form rather quickly, making the close correspondence between 
the results likely in practice. 



C. Idealized collision scenario 



The final test situation that we consider, potentially the most challenging for a description 
such as SHE, is the early stage of a heavy-ion collision, where interpenetration of the opposing 
nuclei has started but no equilibration has yet occurred. To represent such a situation, Welke 
et al. took a momentum distribution with two Fermi spheres, each representing the ground- 
state of saturated matter 

/(r, P ) = ^^%t 0) -|pTPo/2|), (22) 

separated by po — 800MeV/c in momentum space. The scenario is sketched in Fig. 11. 
To account for relativistic effects in this colliding nuclear matter scenario, one would rather 
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FIG. 7. Optical potential for an isotropic momentum distribution at p = 2po, in the local cm., for 
different indicated momenta p. 
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FIG. 8. Optical potential for an anisotropic momentum distribution with e = 0.5 at p = 2po, in 
the local cm., for different indicated momenta p. 



consider Fermi ellipsoids than spheres [13-15]. Such a Lorentz covariant treatment goes 
beyond the goals of the present paper. The optical potentials U (cf. Fig. 12) associated 
with it are plotted as functions of the polar angle & (defined in Fig. 11) for different 
momenta. The vector momenta p are taken with reference to the center of one sphere, with 
the angle relative to the symmetry axis of the system. The mean fields vary not only in 
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FIG. 9. Optical potential for an anisotropic momentum distribution with s = 1.0 at p = 2po, in 
the local cm., for different indicated momenta p. 
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FIG. 10. Optical potential for an anisotropic momentum distribution with e = 1.5 at p = 2po, in 
the local cm., for different indicated momenta p. 



overall offset, but also in dependence on angle. As a point of view additional to that in 
Welke et al. [9], presented above, we further consider the mean- field potential at different 
p relative to the cm., as a function of angle i? about the cm. (defined in Fig. 13). The 
SHE, WPKDG and GBD potentials are compared in this fashion in Fig. 14. Significant 
deviations can be seen between the potentials at lower momenta p, but lesser at the highest 
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FIG. 11. Situation within the local rest frame of two Fermi spheres separated by po in momentum 
space. The center of the left sphere serves as point of origin in the first examination of U for that 
situation. 





FIG. 12. Optical potential for the situation represented in Fig. 11, for different momenta relative 
to the center of one of the Fermi spheres, as a function of angle about that center, for the different 
parametrizations of the potential. 
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FIG. 13. Situation within the local rest frame of two Fermi spheres separated by po in momentum 
space. The center of mass serves as point of origin in the second examination of U for that situation. 

momentum towards the region where the nucleons are. By construction, GBD potential 
is isotropic in this representation. The discrepancies at low momenta in Fig. 14 show 
limitations of the SHE approach. While we chose another phenomenological approach as 
reference here, rather than a microscopic theory for which sparse non-equilibrium results 
exist, we expect a difficulty for SHE, with just a couple of separable terms, to be universal 
when momentum distributions change abruptly with momentum. Fortunately, in heavy-ion 
collisions the distributions quickly evolve to a smooth form. 



IV. DISCUSSION AND CONCLUSIONS 
A. Application to transport simulations 

Beyond the scope of this article is the actual application of the anisotropic mean field 
parametrization to reaction simulations, namely BUU calculations. However, in order to 
stress the practical relevance of our work, we will provide a brief prescription how to employ 
our results in reactions and name the benefits one should gain then. 

The foundation of BUU methodology is the Boltzmann transport equation [12, 21] for the 
single-particle distribution function /(p, r, t). In order to determine the state of the system 
one needs to solve the BUU equation for / in every time step t: 

?L + ^1^1 _ = j (/) . (23) 

dt dp dr dr dp °° 

The collision integral on the r.h.s. of (23) accounts for the time evolution of the single- 
particle phase-space density / evoked by two-body collisions. The mean field U enters the 
Boltzmann equation on the l.h.s. which takes the motion of particles in the field into account. 
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FIG. 14. Optical potential for the situation represented in Fig. 13, for different momenta relative to 
the local center of mass, as a function of angle about the center, for the different parametrizations 
of the potential. 



In order to find the single-particle energies e one has to take the functional derivative of the 
system's net energy E with respect to the phase-space density /: 

5E 

f_ 6V_ 

U . (24) 
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If the p-dependent part of the mean field is parametrized in terms of an integral 

f/(r,p)oc / rfV /( ^ l2 , (25) 

as in the model by Welke et al, one has to determine a different 3-dimensional integral for 
every momentum p and position r. The latter is either computationally costly or gives rise 
to large errors for a low number of Monte Carlo sample points in phase space, corresponding 
to so-called test particles [21]. In the framework of our SHE model one would have to 
evaluate the integrals in the isotropic term and in T zz once per time step for a given spatial 
location. This is a significant advantage of our model over the WPKDG parametrization. 
Since the primary anisotropy occurring during much of the reaction is that of a broadened 
/ in the beam as compared to transverse direction, our model - where an enhancement of 
T zz , as compared to T xx or T yy , is intrinsically built-in - naturally qualifies for describing 
such systems in an effective way. 
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B. Summary 

We developed a parametrization for the nucleonic mean field U, which explicitly exhibits 
an anisotropic behavior for anisotropic phase-space densities /, and that can significantly re- 
duce computational effort and statistical noise in transport simulations, compared to current 
practice. On that account, we made the nuclear energy functional separable in momentum 
space, with different terms corresponding to different multipolarity in spherical angle. As 
reference and guideline during the process of setting up our model, the parametrization by 
Welke et al. [9] was used. We step by step evolved the elements of our model by breaking 
up the original - and in comparison to our parametrization very costly - convolution within 
the potential energy density, substituting it with separable scalar and tensorial terms, both 
symmetric in the single-particle properties, and by taking the functional derivative of V with 
respect to / to obtain U in a simple form. In addition to the Skyrme-type parameters in 
the p-dependent part, our model comprises two exclusive parameters, one of them - A iso - 
mainly relevant in isotropic scenarios, the other one - A an ; so - important for characterizing 
the anisotropy of the mean field. The latter phenomenological parameters represent range 
and are physically expected to be similar. In the framework of our SHE model we can 
reasonably well describe cold nuclear matter properties, obtain excellent results for equili- 
brated scenarios and find a good practical agreement for early stages of a colliding system. 
The latter situation has been dealt with in the literature [9]. We showed that our ansatz is 
generally applicable even when momentum population for the system changes abruptly. 

To conclude, we presented a simple way of describing the potential energy density V and 
associated anisotropic momentum dependent mean field U in heavy-ion collisions. Flexibility 
in coping with interactions for anisotropic momentum distributions is important since local 
momentum anisotropies persist until late in the dynamics of heavy-ion collisions. Explo- 
ration of consequences of the anisotropies, however, should not be stacked against inabilities 
to carry out collision calculations, poor statics or noise in the calculations. With the pro- 
posed approach, we believe, we resolve these issues. 
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Appendix: Analytical expressions for ground-state quantities 

When dealing with Fermi spheres in cold matter scenarios one can find analytical results 
for the integrals over a Yukawa interaction kernel, appearing in (5), (6) and (11), for instance. 
Thus, one can circumvent many-dimensional numerical integrations. This obviously does 
not help in carrying out collision simulations, as the distributions quickly depart from the 
Fermi spheres. Particularly useful are the following identities: 
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